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PERFORMANCE  LIMITS  OF  A JOSEPHSON  JUNCTION  MIXER 


J.  H.  Claassen  and  P.  L.  Richards 

Department  of  Physics,  University  of  California 
Berkeley,  California  94720 

ABSTRACT 

We  present  the  results  of  extensive  analog  computer  simulations  of 
a Josephson  junction  used  as  a mixer  with  an  external  local  oscillator. 
The  Resistively  Shunted  Junction  model  was  used  throughout,  and  the 
effects  of  intrinsic  junction  noise  were  included.  When  the  source 
impedance  is  much  greater  than  the  junction  resistance  R our  simulations 
permit  predictions  of  conversion  efficiency,  noise  temperature,  and 
saturation  level  to  be  made  for  a wide  range  of  experimental  parameters. 
The  possibility  of  harmonic  mixing  has  also  been  considered.  With  a 
source  resistance  comparable  to  R (i.e.,  conventional  "matching")  the 
behavior  of  the  system  is  too  complicated  to  permit  such  general  predic- 
tions of  performance.  From  the  results  of  simulations  for  a particular 
set  of  parameters  it  appears  that  the  best  noise  temperature  will  usually 
be  achieved  for  a source  Impedance  somewhat  greater  than  R.  An  upper 
limit  for  the  mixer  noise  temperature  is  'v  40  T over  a broad  range  of 
junction  parameters,  where  T is  the  effective  temperature  of  the  junction 
The  Conversion  efficiency  under  these  circumstances  should  be  comparable 
to  what  is  potentially  available  from  Schottky  diode  mixers.  Our  simula- 
tions show  that  with  stronger  microwave  coupling  it  is  possible  to  obtain 
conversion  gain  - i.e.,  a conversion  efficiency  exceeding  unity.  This, 
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however,  will  probably  be  at  the  expense  of  a higher  mixer  noise  tempera- 
ture. Saturation  can  be  important  even  for  source  temperatures  'V'  300  K. 

To  avoid  saturation  it  is  necessary  either  to  restrict  the  coupling 
bandwidth  or  use  an  array  of  junctions. 

ij 
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I.  INTRODUCTION 


An  externally  supplied  rf  current  has  a strong  effect  on  the  static 
current-voltage  (I-V)  characteristic  of  a Josephson  junction.  This 
dependence  Is  sufficient  to  make  the  Josephson  junction  a candidate  as  a 
mixer  element  for  use  with  an  external  local  oscillator  (LO) . The  sum 
of  a small  signal  current  and  a relatively  larger  LO  current  is  equiva- 
lent to  an  amplitude  modulated  LO  current  If  the  signal  Is  small  enough 
that  the  phase  modulation  can  be  neglected.  When  the  junction  Is  biased 
in  such  a way  that  changes  in  its  I-V  curve  are  observed,  a response 
will  be  measured  at  the  modulation  frequency,  which  appears  at  the  frequency 
difference  (IF)  between  the  signal  and  the  LO.  If  the  signal  Is  suffi- 
ciently small,  the  output  amplitude  should  be  proportional  to  signal 
amplitude.  This  type  of  mixing  In  a Josephson  junction  was  first  observed 
by  Grimes  and  Shapiro^  In  1968. 

Because  of  their  high  response  speed  and  low  operating  temperature, 
Josephson  junction  mixers  show  particular  promise  In  the  mm  wave  region 
as  front  ends  for  low  noise  receivers.  The  most  often  quoted  figure  of 
merit  for  a receiver  Is  Its  single  sideband  system  noise  temperature: 


T 

sys 


L(T^  + Tj^/n). 


(1) 


L Is  the  loss  between  the  system  Input  and  the  mixer,  which  can  in  prin- 
ciple be  made  close  to  unity,  T^^  Is  the  noise  temperature  of  the  IF 

(difference  frequency)  amplifier  following  the  mixer,  and  T Is  the  noise 

n 

contribution  of  the  mixer  referred  to  the  mixer  Input.  The  conversion 
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efficiency  n is  defined  as  the  ratio  of  power  delivered  to  the  IF  ampli- 
fier, to  the  incident  signal  power.  Evidently  is  the  smallest  possible 
system  noise  temperature  that  can  be  achieved  by  reducing  the  input  losses 
and  the  noise  of  the  IF  amplifier,  and  as  such  represents 
a limit  on  performance  set  by  the  mixer.  Other  mixer  parameters  that 
can  be  important  are  the  saturation  level  and  the  Instantaneous  bandwidth. 
The  latter  is  defined  as  the  range  of  frequencies  around  the  LO  that  can 
be  down-converted. 

In  this  paper  we  will  develop  predictions  of  the  Important  figures 
of  merit  for  Josephson  junction  mixers,  and  determine  the  junction  param- 
eters, bias  conditions,  coupling  circuits,  etc.,  that  optimize  their 
performance.  The  mathematical  model  of  a junction  (the  Resistlvely 
Shunted  Junction  or  RSJ  model)  that  was  used  in  our  calculations  is 
presented  in  section  II  and  Important  scaling  parameters  are  identified. 
Section  III  summarizes  computer  simulations  of  junction  I-V  curves  under 
conditions  of  weak  rf  coupling.  In  section  IV  an  expression  for  conver- 
sion efficiency  is  obtained,  both  for  fundamental  and  harmonic  mixing. 

It  is  found  that  n can  exceed  unity  under  some  circumstances,  and  for 
ideal  junctions  should  have  useful  values  well  into  the  submillimeter 
region.  We  then  discuss  in  section  V the  noise  model  that  is  appropriate, 
and  give  the  results  of  computer  simulations  for  output  noise  and  the 
resulting  mixer  noise  temperature  of  weakly  coupled  junctions.  A simple 

expression,  T < 40  T,  is  found  to  apply  for  a broad  range  of  parameters, 
m 

Saturation  effects  are  estimated  in  section  VI.  We  find  that  such  effects 
limit  the  bandwidth  of  Josephson  junction  mixers  which  are  coupled  to 
room  temperature  sources.  The  analysis  given  in  this  paper  builds  upon 
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concepts  developed  In  previous  publications  of  the  authors  along  with 

2 3 4 

their  co-worker  Y.  Taur.  * * An  attempt  has  been  made,  however,  to 

make  this  discussion  self-contained. 

When  an  rf  source  Is  strongly  coupled  to  a junction,  the  computer 

predictions  of  n and  T depend  In  a complicated  way  on  Junction  par2uneters. 

m 

In  section  VII  we  summarize  the  trends  of  these  dependences  for  specific 
choices  of  parameters.  It  becomes  clear  that  relatively  slight  Improve- 
ments are  achieved  by  strong  rf  coupling,  and  then  only  If  rather  stringent 
requirements  on  junction  parameters  can  be  met.  In  section  VIII  we 
explore  the  possible  use  of  arrays  of  junctions  as  mixers.  Finally, 
we  conclude  In  section  IX  with  a discussion  of  the  properties  of  real 
junctions  that  are  currently  available  and  that  can  be  expected  to 
become  available.  The  effects  on  our  predictions  of  self-heating  and 
deviations  from  the  RSJ  model  are  estimated. 

II . ELECTRICAL  MODEL 

The  Reslstively  Shunted  Junction  (RSJ)  model  of  a Josephson  junction 
that  we  have  used  In  our  calculations.  Is  shown  In  Fig.  1(a),  The  Josephson 
tunneling  current  path  is  represented  as  being  in  parallel  with  a linear 
resistance  R,  and  a capacitance  C.  The  pair  tunneling  current  Is  given 
by  the  well  known  expression,^ 

V(t’)dt' j , 


I ■ I sin 
P c 


(2) 
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where  V is  the  Instantaneous  voltage  across  the  junction.  A theoretical 
expression  connects  R with  the  critical  current  for  at  least  two  kinds 
of  Josephson  junctions, 

RI  - TiA^/4ekT  . (3) 

c c 

Here  A is  the  energy  gap  of  the  superconductor  from  which  the  junction 

is  made,  and  is  its  transition  temperature.  The  data  from  available 

junctions  of  many  types  suggest  that  (3)  gives  an  approximate  universal 

upper  limit  for  the  product  of  R and  that  can  be  obtained.  On  the 

other  hand,  the  shunt  capacitance  C can  vary  considerably,  from  negli- 

0 

gibly  small  in  mlcrobrldges  to  quite  large  in  sandwich-type  tunnel 
9 

junctions. 

In  the  theoretical  analysis  of  the  RSJ  model  it  is  helpful  to  use 
the  following  reduced  variables: 

Current:  i » I/l 

c 

Voltage:  v = 

Frequency:  = hv/2eRI^ 

Time:  t = t(2eRI  )/h 
c 

Capacitance:  6 = (2eA)R^I  C 

c c 

Impedance:  z = Z/R 
Noise:  T “ 2ekTA>I 

c 

In  the  last  expression,  T is  the  physical  junction  temperature  only  if 
the  noise  currents  in  the  junction  arise  primarily  from  classical  Johnson 
noise  in  the  shunt  resistor.  If  other  sources  of  noise  are  Important 
it  must  be  considered  to  be  an  effective  noise  temperature. 
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Practical  constraints  on  the  Impedance  scale  are  Imposed  by  external 
coupling  considerations.  The  characteristic  Impedance  levels  at  both 
the  rf  and  IF  frequencies  must  In  practice  be  of  the  order  of  ten  to  a 
few  hundred  ohms.  For  most  experimental  conditions,  the  junction  resis- 
tance R Is  considerably  less  than  the  magnitude  of  any  external  Impedances 
connected  to  the  junction,  so  that  It  Is  appropriate  to  consider  all 
externally  applied  signals  as  current  sources.  The  case  of  strongly  , 

coupled  junctions  will  be  discussed  In  section  VII.  If  1^  represents 
the  sum  of  these  externally  applied  currents,  the  equations  which 
represent  the  RSJ  model  of  Fig.  1(a)  are; 


V 


d(i> 

dT  ’ 


(4) 


1 

X 


+ 1 

n 


sin  <()  + V + 3 


dv 
c dx 


(5) 


Here  <()  Is  the  difference  between  the  phases  of  the  superconducting  order 

parameter  on  either  side  of  the  junction.  The  noise  current  1^  Is  usually 

2 

assumed  to  have  a "white"  spectral  density,  with  1^  = — FAS).  Here  A(2  Is 

n TT 

the  bandwidth  over  which  the  noise  Is  observed.  However,  for  very  low 
temperatures  a different  noise  spectrum  should  be  assumed,  as  will  be 
discussed  In  a later  section, 

A considerable  fraction  of  the  body  of  literature  on  Josephson 
junctions  consists  of  solutions  to  Eqs.  (4)  and  (5)  under  various  condi- 
tions. A comprehensive  summary  of  the  results  of  these  calculations  Is 
given  In  Ref.  10.  It  should  be  noted  that  only  for  certain  special  cases 
are  analytic  solutions  possible.  The  analog  computer  has  been  a powerful 


J 


tool  in  analyzing  the  Josephson  equations  since  reasonably  accurate 
solutions  can  be  rapidly  generated  as  various  parameters  are  altered. 

The  theoretical  results  in  this  paper  were  all  obtained  from  one  or 

11  12 

another  of  the  junction  simulator  schemes  that  have  been  published.  ’ 


III.  RESULTS  OF  MODEL  CALCULATIONS 

For  > 0.8  the  static  I-V  curve  with  no  rf  current  applied  becomes 

3 

hysteretic  in  the  region  near  zero  voltage  when  F = 0.  The  current 
spread  of  the  hysteresis  loop  increases  monotonically  with  g^.  In  the 
presence  of  significant  rf  current,  a wide  variety  of  effects  occur, 
including  the  possibility  of  hysteretic  I-V  curves.  Hysteretic  junctions 

3 

have  not  generally  been  used  as  mixers,  since  the  appropriate  bias  point  (to  be 
discussed  in  more  detail  in  a subsequent  section)  is  not  readily  acces- 
sible. Sandwich-type  tunnel  junctions  have  not  been  used  as  mixing 
elements  since  these  junctions  have  historically  had  very  large  values 
of  g^,  and  are  strongly  hysteretic.  The  large  capacitance  also  gives  a 
small  rf  impedance  which  leads  to  coupling  difficulties. 

We  shall  consider  the  case  of  a junction  with  g^  = 0 which  is 
driven  by  dc  and  rf  current  sources  as  a first  approximation  to  the 
conditions  under  which  a Josephson  junction  is  actually  used  as  a mixer. 
In  reduced  units  we  write  i * i,  + i ..  sin  fix.  Figs.  2(a)  and  (b) 
show  representative  I-V  curves  with  (2  < 1 for  various  values  of  rf 
current.  Evidently  when  F * 0 we  can  characterize  these  curves  by  the 
dependence  of  the  various  step  heights  on  rf  current.  Unlike  the  case 


> 1,  no  analytical  expressions  for  the  step  heights  exist  when  n < 1, 
except  in  a perturbation  theory  approximation.^^  Various  authors^^’^^’^^ 
have  published  computer  calculations  of  the  step  heights  for  values  of 

n < 1. 

In  this  paper  we  limit  our  attention  to  the  range  0.1  < < 1.0. 

For  normalized  frequencies  greater  than  one  there  is  a rapid  decline  in 

mixer  performance  (q  T On  the  other  hand  we  shall  find 

m 

that  performance  is  not  especially  improved  as  is  reduced  below  0.1. 

The  normalized  frequency  range  we  consider  can  correspond  to  quite  high 
real  frequencies.  For  Nb  Junctions  whose  RI^  product  approaches  the 
theoretical  limit  (3),  for  example,  = 1 corresponds  to  v **  1000  GHz. 
From  the  point  of  view  of  mixing, the  dependence  of  the  0th  step  i 

o 

on  rf  current  is  of  particular  interest.  We  note  from  Fig.  2(d)  that 
there  is  a considerable  range  over  which  the  dc  current  depencs  linearly 
on  the  rf  current.  An  ac  current  proportional  to  di^/di^^  appears  in 
the  output  when  the  LO  current  is  amplitude  modulated  at  S),®*  2(d) 

shows  that  a LO  bias  current  which  reduces  the  0th  step  by  50-80%  will 
produce  the  maximum  IF  output  current. 

In  order  to  consider  the  case  of  harmonic  mixing  (where  the  detected 
signal  is  near  N^q»  with  N an  integer),  we  have  generalized  the  slope 
parameter  dl^/dl^^  by  considering  a combination  of  currents, 

Irf  = ij^Q  sin  + 1^  sin  (N(h-  + 0),  (6) 

where  the  signal  current  1 « i . With  the  simulator  we  found  that  in 


the  range  0 < v < fi,  the  effect  of  i^^  is  to  change  the  current  by  an  amount 
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Ai  = - i.  cos  (6  + 6„)  + 0(1*,).  (7) 

O 1 N N 1 

(Similar  changes  occur  In  other  regions  of  the  I-V  curve,  but  were  not 

Investigated.)  Here  Is  a generalized  slope  parameter,  with  = 

dl  /dl, and  6 =0.  With  the  simulator  we  find  that  all  the  S.,  are 

o LO  1 N 

nearly  constant  over  the  range  of  LO  currents  which  reduce  the  zero 

voltage  current  by  a factor  of  order  2,  and  6^  = 90°,  6^  = 180°,  etc. 

In  Fig.  3 we  show  the  dependence  of  the  first  few  on  normalized  fre- 
quency as  determined  by  the  simulator,  for  a value  of  LO  current  that 

reduces  the  critical  current  to  0.4  I . 

c 

A rather  different  mode  of  harmonic  mixing  for  even  values  of  N 
2 

has  been  observed.  When  the  LO  current  Is  adjusted  to  completely 
suppress  the  zero  voltage  current,  efficient  mixing  occurs  at  zero  dc 
bias.  The  analysis  of  this  mode  of  operation  requires  a different 
approach  than  that  presented  here,  so  Is  outside  the  scope  of  the 
present  paper. 

We  can  now  consider  the  combination  of  a LO  current  and  a signal 

current  which  differs  In  frequency  by  a small  amount  (2  from  or 

IF  LO 

one  of  Its  harmonics,  by  setting  9 = Ot-T  In  Eq.  (6),  Then  according 

Xi* 

to  (7)  there  Is  an  output  current  at  the  difference  frequency  0 , 

IF  IF 

with  an  amplitude  1,S„. 

1 N 

Since  the  quantity  of  Interest  Is  the  available  output  power , we 
must  consider  the  output  Impedance  as  well  as  the  output  current.  The 
output  Impedance  Is  the  differential  Impedance  at  the  bias  point. 

As  Is  Indicated  In  Fig.  2(b),  the  differential  resistance  between  steps 
Is  practically  Infinite  for  (2  < 0.6  when  there  Is  no  noise.  When  noise 
Is  added,  the  corners  of  the  steps  are  rounded  and  a maximum  In  occurs 
midway  between  steps  for  typical  values  of  F.  With  the  simulator  we  find 
that  an  empirical  expression  for  r^  midway  between  the  zero’th  and  first 
step  Is: 

tp  - p(r  + r^)“** . (8) 

The  dependence  of  P and  F^  on  normalized  frequency  Is  shown  In  Fig.  4. 
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When  r < a somewhat  higher  value  of  than  given  by  (8)  can  be  found 
near  a step.  For  achievable  values  of  F this  is  only  likely  to  occur 
when  ^ ^ 1,  and  throughout  this  paper  we  assume  a bias  at  the  midpoint. 

As  can  be  seen  from  Fig.  2,  the  static  I-V  curves  for  a progression  of 
LO  levels  are  parallel.  This  means  that  the  output  impedance  will  be 
insensitive  to  LO  drifts  for  a constant-voltage  biasing  scheme. 

It  should  be  mentioned  that  in  practice  the  differential  resistance 
is  a parameter  that  seems  to  be  quite  sensitive  to  deviations  from  our 
model  caused  by  effects  such  as  shunt  capacitance,  spurious  external 
resonances,  and  various  kinds  of  structure  in  the  I-V  curve  that  cannot 
be  explained  by  the  RSJ  model.  The  latter  Includes  energy  gap  structure, 
"excess  current",  etc.  In  comparing  experimental  results  with  theory 
it  is  often  appropriate  to  treat  r^  as  a measured  parameter.  One  should 
expect,  however,  that  the  dependence  implied  by  Eq.  (8)  will  be  obeyed 
qualitatively,  i.e.  r^  is  reduced  when  I^  is  reduced  or  T is  Increased,  etc. 

The  equivalent  output  circuit  of  Fig.  1(c)  certainly  applies  in  the 
limit  of  sufficiently  low  difference  frequencies.  The  upper  frequency 
limit  for  which  it  is  applicable  depends  on  the  speed  with  which  the 
phase  evolution  of  the  junction  approaches  the  steady-state  solution 
after  a sudden  change  in  When  the  bias  point  is  roughly  midway 

between  steps,  this  should  occur  after  only  a few  LO  cycles. We  have 
directly  tested  the  frequency  dependence  of  the  output  current  and 
impedance  on  a junction  simulator  for  the  case  F = 0.005,  (2  = 0.5,  and 
find  that  the  low  frequency  model  remains  quite  accurate  up  to  (2  /(2  =0.1. 

Even  when  " 0.3,  the  deviation  from  the  low  frequency  model  is 

less  than  '''  20%. 


1 
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IV.  CONVERSION  EFFICIENCY 


An  Important  figure  of  merit  for  mixers  is  the  conversion  efficiency, 
or  "gain",  defined  as  the  ratio  of  power  delivered  to  the  IF  amplifier  to 
the  available  rf  signal  power.  An  expression  that  is  applicable  to  any* 


nonlinear  system  is 


P 


Here  P is  defined  as  the  available  LO  power, 

lAJ 

and  Cjp  is  the  coupling  efficiency  to  the  IF  amplifier  input  impedance  R^, 
given  by  C^^  = 4 ^ what  follows  we  will  assume  that  the 

IF  amplifier  is  matched  to  the  junction  (R^^  = R^^)  so  that  C^.^  = 1. 

It  may  be  argued  that  connecting  an  IF  amplifier  of  finite  input 
Impedance  to  the  junction  violates  our  assumption  of  external  current 
sources.  It  can  be  demonstrated,  however, that  the  dc  I-V  curve  of  a 
junction  depends  only  on  the  external  Impedance  at  frequencies  of  the 


order  or  greater  than 


the  Josephson  oscillation  frequency . Any 


reasonable  IF  amplifier  coupling  would  be  effective  only  at  much  lower 
frequencies;  for  our  simulations  we  model  the  high  frequency  decoupling 
by  a series  Inductance  as  is  shown  in  Fig.  1(b). 

The  series  resonance  shown  in  Fig.  1(b)  is  the  simplest  representa- 
tion of  the  rf  coupling  circuit  that  would  apply  to  a real  experiment. 


For  a nonlinear  reactance'*’  Eq.  (9)  would  n«c  correctly  predict  the  con- 
version efficiency,  which  is  proportional  to  the  IF  frequency.  The 
expression  given  is  a lower  limit  to  r)  for  an  arbitrary  device. 
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This  reflects  the  fact  that  useable  Josephson  junctions  have  a shunt 

resistance  that  Is  small  compared  to  typical  antenna  Impedances,  which 

19 

are  of  the  order  of  several  hundred  ohms,  except  at  particular  resonant 
frequencies.  Our  simulations  show  that  with  resonant  coupling  of  the  LO, 
the  junction  response  can  be  fairly  well  fit  by  the  expression. 


Here  vdilch  Is  shown  as  a function  of  normalized  frequency  In  Ref.  3, 

Is  a parameter  that  plays  the  role  of  an  Input  resistance  at  the 

rf  frequency.  The  expression  for  conversion  efficiency  then  becomes 


^ "rf'' 


Optimization  of  n by  varying  R Is  complicated  by  the  fact  that  when  R 

o o 

Is  of  the  order  of  R the  differential  resistance  R^  can  be  strongly 
reduced  from  the  value  given  by  Eq.  (8).  This  effect  Is  particularly 
pronounced  for  small  Thus  the  value  of  R for  which  n Is  maximum  Is 
greater  than  This  somewhat  disturbing  result  occurs  because  a 

Josephson  junction  does  not  behave  as  a simple  two-port  device.  Its 
properties  are  affected  by  the  Impedance  It  sees  not  only  at  the  Input 
and  output  frequencies,  but  at  many  other  frequencies  as  well.  When  the 


noise  performance  Is  considered,  the  optimum  value  of  R can  be  still 

S 

larger.  Hence  the  parameter  R^^  Is  of  little  significance  In  establishing 
the  optimum  source  resistance.  We  put  off  a detailed  discussion  of  these 
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effects  of  strong  coupling  to  a later  section,  and  consider  the  current- 
driven  limit  R » R.  Then  Eq.  (11)  becomes 


2 

n = S?  — , r„  » 1. 

1 's  ^ 


In  the  expression  (12)  for  conversion  efficiency  there  Is  an  Implicit 
dependence  on  the  noise  parameter  F via  the  factor  r^^.  We  have  demons- 
trated with  the  simulator  that  Is  Independent  of  F.  Although  the 

values  of  S for  N > 1 were  obtained  with  F ■ 0 It  Is  plausible  that  they 
N 

too  are  Independent  of  F.  The  dependence  of  conversion  efficiency  on 
Junction  parameters  for  a given  normalized  frequency  and  rf  coupling, 

-jj  u 

n '''  F”^  'v<  (I  /T)^,  suggests  that  I should  be  large.  Since  the  RI 
c c c 

product  Is  roughly  constant,  this  conflicts  with  a practical  requirement 

of  microwave  matching,  that  R should  not  be  < 10  ohm. 

For  the  case  of  harmonic  mixing  ~ ^ ^ conversion 

efficiency  Is  given  by  Eq.  (12)  with  S,  replaced  by  S„.  It  Is  some 

1 N 

Interest  to  compare  the  expected  mixing  efficiency  for  various  orders 

keeping  the  signal  frequency  fixed,  l.e.  using  a LO  frequency  near  (2  ^ , 

slg 

^gj^g/3.  etc.  In  Fig.  5 we  plot  which  Is  proportional  to 

the  conversion  efficiency,  as  a function  of  . and  N.  Note  that  S„ 

slg  N 

and  r„  are  evaluated  at  ^ /N.  We  have  assumed  F = 0.005,  as  would 
D slg 

result  for  Instance  If  T * 4.2  K,  I “35  pA-  The  precipitous  decline 

c 

In  conversion  efficiency  for  the  fundamental  mixing  case  above  J2  'V'  0.6 
Is  somewhat  exaggerated  due  to  the  assumption  of  a bias  point  midway 
between  the  zeroth  and  first  steps.  In  fact,  a bias  point  with  higher  r 


can  usually  be  found  closer  to  one  of  the  steps.  This  becomes  Increasingly 
true  for  fi  > 1. 


In  a recent  experiment  by  D'yakov  e£  on  point  contact  Josephson 

Junctions  poorly  coupled  to  the  microwave  source  (and  thus  satisfying 
the  criterion  of  being  current-driven)  the  observed  fall  off  with  harmonic 
order  was  much  slower  than  indicated  in  Fig.  5.  We  find  that  most  of  the 
discrepancy  comes  from  deviations  in  the  experimental  values  of  r^  from 
the  predictions  of  the  resistive  model.  The  observed  values  of  the  ratio 
S3/S1  were  generally  somewhat  greater  than  our  prediction,  in  some  cases 
by  as  much  as  35%. 

The  fact  that  Josephson  effect  mixers  require  very  little  LO  power 
means  that  it  should  usually  be  possible  to  derive  adequate  power  for 
fundamental  mixing  from  conventional  sources  using  diode  harmonic  genera- 
tors. It  therefore  seems  that  it  would  be  worthwhile  to  make  use  of  the 

relatively  efficient  harmonic  mixing  capabilities  of  the  Josephson  device 

21  22 

only  for  special  applications  such  as  frequency  comparison,  ’ etc. 


V.  MIXER  NOISE 


The  second  important  figure  of  merit  for  a mixer  is  its  single  side- 
band noise  temperature  T^.  We  obtain  this  quantity  by  first  computing  an 

output  noise  temperature  T which  is  related  to  T by  T * T /n.  To 

out  m m out 

predict  the  output  noise  we  need  a model  for  the  equivalent  driving  noise 

current  1^  of  Eq.  (5).  We  consider  an  expression  that  was  first  given  by 
23 

Scalaplno  for  tunnel  junctions: 

l“(v)  - + ^)coth  + I^(V  - (13) 
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Here  V is  the  frequency  at  which  the  noise  is  observed,  B is  the  band- 
width, and  Iq(v)  is  the  quaslpartlcle  current  in  the  junction  at  a bias 
voltage  V.  When  is  proportional  to  V,  as  is  assumed  in  the  RSJ  model, 
the  above  expression  reduces  to 

I^(v)  = ^ (14) 

where  E is  the  largest  of  eV,  hv,  or  2kT.  We  will  identify  the  regimes 
corresponding  to  the  first,  second,  or  third  term  being  dominant  as  shot 
noise,  photon  noise,  or  thermal  noise,  respectively.  Almost  all  discus- 
sions in  the  literature  of  noise  effects  in  the  RSJ  model  have  assumed 
the  thermal  noise  limit.  The  parameter  F evaluated  at  the  physical 
temperature  T characterizes  the  driving  noise  amplitude  for  this  case. 

In  general,  the  driving  noise  is  processed  by  various  nonlinear 
interactions  within  the  junction  to  produce  a low  frequency  voltage  noise 
that  is  greater  than  RjjI^(O) . Likharev  and  Semenov  calculated  the 

voltage  noise  response  to  a driving  current  noise  I for  the  RSJ  model 

n 

without  rf  bias: 

^n^°^  “ I^(2eV/h)/2i''j  , (15) 

where  i is  the  normalized  bias  current.  This  expression  is  valid  as  long 
as  noise  rounding  is  not  severe. 

Experimental  results  for  the  voltage  noise  across  point  contacts 

25 

without  rf  bias  are  not  consistent.  Ranter  and  Vernon  observed  noise 
at  'V'  150  KHz,  and  the  data  presented  are  considerably  in  excess  of  Eq.  (15), 
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especially  at  low  bias  voltage.  They  note,  however,  that  some  junctions 
showed  a dependence  on  voltage  that  is  qualitatively  similar  to  Eq.  (15). 
Insufficient  Information  is  given  about  these  junctions  to  make  quantita- 
tive comparisons.  Noise  measurements  in  our  laboratory  on  oxidized  Nb 

4 * 

point  contacts  at  50  MHz  are  in  qualitative  agreement  with  Eq.  (15), 

but  quantitatively  in  excess  by  approximately  a factor  2. 

Measurements  of  the  llnewldth  of  the  Josephson  oscillation,  which 

should  be  proportional  to  the  low  frequency  voltage  noise,  have  been 
26  27 

made  both  directly  and  indirectly  in  point  contacts.  At  high  bias 


voltages  these  results  scale  as  would  be  expected  from  Eq.  (15)  in  the 
appropriate  shot  noise  limit,  but  with  a magnitude  a factor  2 or  so  too 


high.  At  low  voltages  the  noise  is  much  higher  than  given  by  Eq.  (15). 

The  results  of  some  of  the  experiments  mentioned  above  were  compared 

28 

with  an  expression  derived  by  Stephen  for  the  voltage  noise  across  a 
tunnel  junction  biased  on  a cavity  mode  step.  Whatever  agreement  was 
observed  would  seem  to  be  fortultuous,  since  the  model  used  for  Stephen's 
calculation  does  not  apply  to  point  contacts.  A certain  amount  of  con- 
fusion has  resulted  from  the  Identification  of  one  of  the  terms  in 
Stephen's  expression  as  "pair  shot  noise".  In  fact,  this  term  arises 
as  a result  of  detection  in  the  Josephson  junction  of  photon  fluctuations 
in  the  cavity,  rather  than  intrinsic  fluctuations  in  the  pair  current. 

As  we  have  defined  it,  shot  noise  in  junctions  for  which  the  RSJ 
model  applies  may  be  considerably  reduced  from  the  value  given  by  Eq.  (13). 

Full  shot  noise  is  expected  when  the  voltage  drop  occurs  across  a region 

29 

much  smaller  than  an  electron  mean  free  path,  as  is  the  case  for  tunnel 


* 

Including  corrections  that  arise  when  noise-rounding  is  severe. 
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* 


junctions.  For  many  low-capacitance  structures,  such  as  the  variable- 
thickness  bridge,  however,  this  may  not  be  the  case. 

There  exist  no  analytic  expressions  In  the  literature  for  the  low 
frequency  noise  output  from  a Josephson  junction  biased  for  optimum  mixer 
operation.  We  have  used  the  analog  simulator  to  determine  the  output 
noise  vdten  it  is  driven  with  an  appropriate  noise  current,  for  a variety 
of  junction  parameters.  Since  at  least  some  of  the  available  experimental 
data  are  consistent  with  calculations  based  on  Eq.  (14)  as  the  correct 
driving  noise,  we  have  used  this  expression  in  the  simulator  calculations 
with  the  philosophy  that  the  results  should  represent  a lower  limit  on 
the  noise  in  a junction. 

There  are  two  extreme  limits  that  can  be  considered:  purely  thermal 

noise  (2kT  » hv  ) and  purely  photon  noise  (hv  » 2kT) . These  corres- 
LiU  iaU 

pond  to  a white  driving  noise  spectrum  and  a spectral  density  proportional 
to  frequency  (blue  noise),  respectively.  Both  extremes  are  likely  to  be 
encountered  in  practical  experiments.  We  note  that  the  crossover  between 
regimes  occurs  at  a frequency  v _ = 42  GHz  when  T = 1 K. 

IjU 

For  the  simulations,  the  white  noise  current  was  obtained  by  ampli- 
fying the  Johnson  noise  in  a resistor  using  a low-noise  amplifier  (PAR  185) . 

This  could  be  converted  to  "blue  noise"  by  passing  it  through  a linear 

u 

network  with  a frequency  response  which  varied  as  n/  v^. 

In  the  mixer  application  the  presence  or  absence  of  shot  noise  is 
largely  irrelevant:  In  the  thermal  noise  limit,  we  have  2kT  » eV,  since 
eV  n.  hv^/4  at  the  optimum  bias  point.  In  the  photon  noise  dominated 
region,  the  shot  noise  contribution  is  Important  only  for  frequencies 
leas  than  shall  show  that  noise  currents  in  this  range  contri- 

bute little  to  the  output  noise. 
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We  first  consider  the  case  of  thermal  noise.  In  this  limit  the 

junction  Is  driven  by  white  noise  with  a spectral  density  given  In 

2 2 

normalized  units  by  1^  • — FAfi.  The  spectral  density  of  the  output 
voltage  noise  from  a junction  simulator  was  determined  with  the  aid 
of  a digital  computer  system,  up  to  a frequency  for  various  values 

of  Important  parameters. 

It  Is  most  convenient  to  consider  the  output  noise  current  I as 

no 

Is  Indicated  In  the  equivalent  circuit  of  Fig.  1(c),  rather  than  the  noise 

voltage.  We  define  a noise  parameter  8^  as  the  ratio  of  to  the 

no 

2 

driving  noise  current  at  low  frequencies: 


8" 


1^  (0)  v*(0;r 

no  _ _n 

1^(0)  V ■ 

n 


(16) 


Fig.  6 shows  some  simulator  predictions  of  8^  as  a function  of  bias 
voltage  for  a junction  biased  with  an  rf  current  source  whose  amplitude 
is  appropriate  for  optimum  mixer  operation. 

To  determine  the  mixer  noise  temperature  we  note  that  the  effective 
output  temperature  is  r^T,  where  T is  the  physical  temperature 

of  the  junction.  Then 


T - T Jr\  8^T  r-/S?  , (17) 

m out  o 1 

where  we  have  used  Eq.  (12)  for  n.  The  only  term  in  this  expression  that 
depends  on  bias  voltage  (if  the  voltage  is  below  the  first  step)  is  8^. 
The  minimum  seen  in  Fig.  6 to  occur  midway  between  the  zeroth  and  first 
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steps  is  a feature  observed  for  all  Q and  for  all  values  of  F small  enough 
to  leave  some  step  structure  visible  in  the  I-V  curve.  Hence  the  bias 
point  for  optimum  mixer  noise  temperature  coincides  with  the  point  of 
maximum  conversion  efficiency  when  Q < 1.  We  have  not  studied  the  case 
of  higher  normalized  frequencies,  where  the  optimum  in  conversion  effi- 
ciency can  lie  closer  to  the  zeroth  step. 

The  simulations  show  that  6^  . is  independent  of  F (for  a given  fi) , 

min 

as  long  as  F is  small  enough  to  leave  r^  > 1.  It  is  also  unchanged  over 

the  range  of  LO  levels  appropriate  for  mixer  operation,  i.e.  those  that 

suppress  the  critical  current  by  50-80%.  Hence  is  a function  only 

of  normalized  frequency.  This  dependence  is  shown  in  Fig.  7.  We  also 

show  which  according  to  Eq.  (17)  should  be  proportional  to  the 

mixer  noise  temperature  in  the  current -biased  limit.  This  factor  is 

practically  Independent  of  normalized  frequency  over  most  of  the  range 

of  Interest  to  us.  Note  that  T does  not  depend  on  F,  in  contrast  to 

m 

the  conversion  efficiency.  This  simplification  unfortunately  does  not 
apply  in  the  case  of  a strongly  coupled  junction  (r  '''  1) , as  we  shall 
show  in  a later  section. 

Eq.  (17)  is  valid  in  the  current-biased  limit,  r^  » 1.  On  the 

basis  of  simulations  with  tg  < 2 to  be  described  in  section  VII,  we  can 

guess  that  (17)  remains  fairly  accurate  with  r as  low  as  4.  Using 

this  value  along  with  the  results  given  in  Fig.  7,  we  find  T^  < 40  T 

for  n in  the  range  0.1  - 0.9.  To  make  a similar  estimate  of  conversion 

efficiency,  we  must  specify  F.  According  to  Fig.  5,  p > 0.25  in  the 

above  range  of  frequencies  if  F < 0.005.  The  latter  condition  would  be 

met,  for  example,  if  I > 35  UA  at  T ■ 4 K. 

c 
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Equation  (17)  suggests  that  can  be  improved  indefinitely  by 
lowering  the  junction  temperature  T.  However,  once  T < hv  /2k.  the 

LiU 

thermal  noise  results  no  longer  apply.  To  investigate  the  low  tempera- 
ture regime,  we  have  driven  the  simulator  with  "blue"  noise  - i.e.  with 
a current  spectrum  given  by  I^(V)  » 2hVB/R.  According  to  Eq.  (13),  this 
corresponds  to  the  limit  T » 0.  In  analogy  to  the  definition  for  , 
a dimensionless  output  noise  parameter  can  be  defined: 


no  n LO 


We  find  again  that  as  long  as  r^^  > 1 this  parameter  is  Independent  of  R. 
Its  dependence  on  normalized  frequency  is  shown  in  Fig.  7.  An  expression 
can  then  be  obtained  for  the  mixer  noise  temperature  in  the  limit  T 0 
which  is  analogous  to  Eq.  (17), 


^ "'■lo 

.2  2k  • 


Since  Y is  of  the  same  order  as  ^ over  much  of  the  range  of  normalized 

min 

frequency  being  studied  here,  we  can  roughly  say  that  in  the  low  tempera- 


ture limit  the  junction  is  at  an  effective  ambient  temperature  hv^^/2k. 

It  is  interesting  to  speculate  on  the  mechanisms  giving  rise  to  an 
equivalent  output  noise  current  near  zero  frequency  greater  than  the 
driving  noise.  One  obvious  source  of  the  extra  noise  is  down-conversion 
of  high  frequency  noise  components  due  to  the  fundamental  and  harmonic 
mixing  processes  we  have  discussed.  This  should  contribute  an  amount 
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00  00 

2 ^ to  , and  2 ^ NS^  to  Y^»  respectively.  The  factor  2 arises 
N-1  ” N-1  " 

from  the  fact  that  two  sidebands  at  each  frequency  are  down-converted 

into  the  IF  bandwidth.  We  find  in  fact  that  the  above  sums  are  cloae 

to  half  of  the  directly  calculated  values  of  B^  - 1 and  respectively. 

30 

We  have  been  advised  that  the  remaining  noise  can  be  accounted  for  by 

processes  analogous  to  mixing  with  the  Josephson  oscillation  In  junctions 
24 

with  no  rf  bias. 


VI . SATURATION 

For  any  mixer  the  conversion  efficiency  will  become  nonlinear  for 
signal  powers  approaching  some  fraction  of  the  LO  power.  In  a Josephson 
effect  mixer  the  LO  level  is  sufficiently  small  that  such  saturation 
effects  could  severely  limit  its  usefulness  in  some  applications.  In 
this  section  we  consider  two  mechanisms  for  saturation. 

One  mechanism  for  nonlinear  response  corresponds  to  output  voltage 
excursions  that  exceed  the  range  of  constant  differential  resistance  on 
the  I-V  curve.  Considering  the  qualitative  I-V  curves  shown  in  Fig.  2, 
it  is  clear  that  when  the  peak-to-peak  voltage  excursions  approach  hv  /2e 

LjU 

31 

the  junction  will  be  strongly  saturated.  More  generally,  we  propose 
that  saturation  occurs  when  the  rms  IF  voltage  exceeds  x*hVj^Q/2e,  where 
X lies  in  the  range  0.1  - 0.2,  depending  on  the  stringency  of  the  criterion 
for  saturation  as  well  as  details  of  the  I-V  curves.  The  corresponding 
signal  power  is 
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P 


sat 


(20) 


In  Ref.  32  data  are  shown  suggesting  that  the  saturation  power  scales  as 
RI^.  According  to  the  above  equation  this  is  true  for  fixed  r^,  and 
rf  coupling.  As  all  of  these  parameters  were  not  recorded,  we  cannot 
make  a detailed  comparison  with  Eq.  (20).  If  we  assume  r 1,  as  was 
typical  of  the  experiments  of  Ref.  32,  the  two  examples  shown  (Fig.  27) 
give  X = 0.23  and  x = 0.33,  respectively,  for  1 dB  gain  compression. 

If  Eq.  (8)  is  used  for  r^^,  with  F » F^,  Eq.  (20)  becomes 


P 

sat 


27rx  fir. 


P^ 


• kTV  > 50  X r kTv. 

LiO  2 


LO 


(21) 


The  numerical  factor  is  valid  for  normalized  frequencies  in  the  range 

0.1  - 1.0. 

A second  saturation  mechanism  is  applicable  to  the  frequently 
encountered  situation  where  the  "signal"  is  in  fact  noise,  e.g.  when  a 
radiometer  is  used  to  observe  the  temperature  T of  a black-body  source 
which  is  coupled  over  a bandwidth  B to  the  mixer.  In  the  limit  of 
current-source  coupling  that  is  being  considered  here,  we  can  say  that 
the  effective  junction  temperature  is  T + T /r  within  the  coupling  band- 
width,  and  T over  the  remaining  frequency  region.  The  differential 
resistance  of  the  junction  will  depend  on  Tg  in  a manner  somewhat  similar 
to  its  dependence  on  T and,  since  the  conversion  efficiency  is  propor- 
tional to  tjj,  the  mixer  will  have  a nonlinear  response  to  changes  in  Tg, 


(22) 


To  estimate  the  magnitude  of  this  effect,  we  use  a generalization  of 
Eq.  (8)  for  the  differential  resistance; 


Here  B is  the  bandwidth  over  which  the  source  is  coupled  to  the  junction, 
and  r = 2ekT  /hi  r . The  dependence  implied  by  Eq.  (22)  was  verified 
for  r = 0 with  the  simulator,  but  not  for  the  general  case  F 0,  F 0. 
The  dependence  of  Q on  normalized  frequency  is  shown  in  Fig.  4.  In  the 
radiometer  mode  one  looks  for  small  changes  in  the  output  power  of 

the  IF  amplifier  corresponding  to  changes  in  the  source  temperature. 
Suppose  in  particular  the  temperature  of  the  source  changes  from  zero 
to  some  value  T . Then  we  have 


AP 


out 


« r (T  )(T  + 2T-)  - r^(0)T  . 
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This  equation  assumes  that  6^,  S*,  and  T^^  are  constant,  a reasonable 

assumption  if  the  IF  ampUEiefris  initially  matched  to  the  junction.  The 
factor  2 in  the  parentheses  results  from  the  fact  that  two  sidebands  are 
down-converted  into  the  IF  bandwidth.  Setting  F^  - 0 for  simplicity, 
we  have 


(T  + 


eT^y 


r*5 


where  c 


(24) 
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As  can  be  seen  in  Fig.  4,  (P/Q)^  is  virtually  Independent  of  normalized 
frequency  and  has  a maximum  value  '^0.45.  It  can  be  shown  that  to  lowest 
order  in  C Eq.  (23)  becomes 


AP  ^ “ 
out 


There  are  two  results  of  the  modulation  of  r^  by  the  source  temperature: 
The  effective  conversion  efficiency  is  reduced  from  what  it  would  be  for 
a coherent  signal,  and  there  is  a square  law  (detection)  contribution  to 
the  IF  response.  The  relative  change  in  conversion  efficiency  is 


l/  B \ ^ 
^VLO/ 
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where  we  have  used  Eq.  (17)  for  T^.  The  requirement  that  this  quantity 
be  small  places  a fairly  stringent  requirement  on  the  allowable  coupling 
bandwidth.  If  the  bandwidth  is  narrow  enough  that  the  conversion  effi- 
ciency is  not  significantly  altered,  the  amount  of  "signal" power  that  will 
result  in  a 1 dB  gain  compression  due  to  the  square  law  term  in  (25)  is 


= 0.2k(f) 


In  practice, care  is  required  to  avoid  confusion  of  these  effects  with 

linear  mixing.  Comparing  with  the  expression  for  output  saturation, 

Eq.  (21),  we  see  that  if  we  choose  x = 0.2,  the  two  mechanisms  are  of 

roughly  equal  Importance  in  determining  a saturation  power.  If  r ■ 4 

o 
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and  T ■ 4 K,  Eq.  (27)  shows  that  a relative  coupling  bandwidth  of  10% 
to  a 300  K source  Is  sufficient  to  compress  the  gain  by  1 dB. 


VII.  STRONGLY  COUPLED  JUNCTIONS 


We  now  consider  the  case  of  strong  rf  coupling,  l.e.  where  the  source 
resistance  R„  Is  of  the  same  order  as  the  Junction  resistance.  When  there 
are  no  reactive  elements  and  the  source  Is  perfectly  coupled  over  the 
complete  bandwidth  from  dc  to  (at  least)  the  gap  frequency,  the  problem 

can  be  cast  Into  the  form  of  a current-biased,  reslstlvely  shunted  junc- 
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tlon.  Not  only  Is  such  a circuit  Impossible  to  realize  In  practice  but. 
In  view  of  the  discussion  In  the  last  section,  saturation  effects  could  be 
expected  to  dominate.  We  have  chosen  to  Investigate  the  series-resonant 
equivalent  circuit  shown  In  Fig.  2(b)  to  represent  the  microwave  coupling 
scheme . 

The  added  complexity  of  such  a circuit  shows  Itself  In  the  appearance 
of  hysteresis  In  the  static  I-V  curve  at  bias  points  appropriate  for  mixer 
operation  when  Rg  Is  of  the  order  of  the  junction  resistance.  An  explana- 
tion for  the  Instability  of  a bias  point  between  the  zero'th  and  first 
steps  Is  given  In  section  VIII.  As  might  be  expected,  the  addition  of 
noise  can  result  In  a "smearing  out"  of  the  hysteretlc  region  to  give 
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a single-valued  dc  I-V  curve.  This  Is  well  known  for  the  case  of  capacl- 
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lively  caused  hysteresis.  In  Fig.  8 we  show  a sequence  of  I-V  curves 
with  various  values  of  the  noise  parameter  F,  obtained  with  a simulation 
of  the  full  circuit  of  Fig.  2(b).  It  Is  noteworthy  that  once  the  noise 
Is  large  enough  to  suppress  the  hysteresis,  the  differential  resistance 
between  the  zero'th  and  first  steps  Is  significantly  less  than  It  would 
be  for  the  current-driven  case.  This  Is  especially  pronounced  at  low 

normalized  frequencies,  where  the  reduction  in  r^  can  be  as  much  as  an 

i 

■ order  of  magnitude  In  the  well  coupled  case. 

I Another  characteristic  of  the  I-V  curves  of  strongly  coupled  junctions 

is  a pronounced  effect  on  the  shape  when  the  LO  frequency  is  detuned  from 

I 

■ the  center  frequency  of  the  resonant  circuit.  An  example  of  this  is  given 

[ in  Fig.  9.  Note  that  by  operating  somewhat  off  resonance  the  differential 

I 

f 

resistance  can  be  increased  above  the  value  at  resonance.  The  complicated 
dependence  of  r^  on  source  resistance,  tuning,  etc.,  makes  it  impossible 
to  arrive  at  a simple  analytical  expression  for  conversion  efficiency  in 
the  strong  coupling  range. 

A second  and  more  important  effect  of  strong  coupling  of  the  LO  is 


r 


4 
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an  increase  in  over  the  value  observed  for  current-source  coupling. 

The  excess  can  be  an  order  of  magnitude  for  typical  experimental  param- 
eters. In  general,  increases  as  r decreases,  becoming  divergent 

as  tp  ■*  <».  Excess  noise  of  this  sort  seems  to  be  a general  feature  of 
hysteretic  I-V  curves  that  have  been  smeared  by  thermal  noise  into  conti- 
nuous ones,  no  matter  what  is  the  mechanism  for  hysteresis.  The  existence 
of  extra  low  frequency  noise  seems  quite  plausible  if  the  smearing  process 
is  regarded  as  one  in  which  the  thermal  noise  randomly  drives  the  junction 
between  its  two  stable  states. 

Since  the  properties  of  strongly  coupled  junctions  depend  on  many 
parameters,  such  as  F and  the  tuning  of  the  resonant  circuit,  that  do  not 
enter  in  the  current-driven  case,  it  is  very  difficult  to  make  general 
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predictions  of  the  performance  to  be  expected.  In  previous  treatments 
we  have  simply  assumed  that  the  noise  with  strong  coupling  is  a factor  2 
or  so  worse  than  with  current  source  coupling.  In  fact,  the  discrepancy 
can  be  much  worse.  We  have  recently  obtained  simulator  predictions  of 
mixing  performance  under  a variety  of  conditions  of  strong  coupling  that 
might  typically  be  used.  These  results, which  are  summarized  in  Fig.  10, 
are  fairly  sensitive  to  details  of  the  model.  Addition  of  a small  shunt 
capacitance  or  a more  realistic  modelling  of  the  rf  circuit,  for  Instance, 
might  have  a substantial  effect.  They  should  not  be  taken  as  exact  numerical 
predictions,  but  as  indications  of  general  trends. 

One  of  the  most  significant  predictions  is  the  rapid  deterioration 

in  T when  I is  Increased,  but  the  RI  product  is  held  fixed.  This 
me  c 

reflects  the  fact  that  0^  , depends  on  F,  in  contrast  to  the  current- 

min 


biased  limit.  In  a mixer  utilizing  a low  resistance  Josephson  structure 
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such  as  a thin  film  weak  link,  the  microwave  circuit  should  be  designed 

for  fairly  weak  coupling  (rg  > 4,  say).  On  the  other  hand,  If  the 

resistance  Is  high  enough  there  may  be  some  advantage  In  reducing  rg  to 

1 or  2.  Note  that  It  Is  almost  never  useful  to  make  r < 1.  Another 

difference  from  the  weakly  coupled  limit  Is  a weaker  dependence  of 

on  the  junction  temperature.  We  have  noted  for  the  former  case  that 

T “ T for  T > hv, -/2k.  For  strongly  coupled  junctions  a temperature 
m LO 

can  generally  be  reached  below  which  no  Improvement  Is  possible,  even  If 

It  Is  still  In  the  thermal  noise  limit. 

We  have  also  considered  the  photon  noise  dominated  limit  for  one  set 

of  parameters  of  a strongly  coupled  junction.  We  find  that  the  results 

for  conversion  efficiency  and  noise  are  the  same  as  If  It  were  driven  by 

white  noise  corresponding  to  F = 2TreV  /I  . If  this  holds  true  generally 

Ijv  c 

we  can  use  the  predictions  of  Fig.  10  In  the  high  frequency  limit  > 2kT 

by  considering  the  junction  to  be  at  an  effective  temperature  hVj^/2k. 

In  the  above  discussion  of  parameter  optimization  we  have  Ignored 
the  variation  of  conversion  efficiency  on  the  grounds  that  It  Is  In  every 
case  large  enough  that  a state-of-the-art  IF  amplifier  will  not  contri- 
bute significantly  to  the  system  noise  temperature.  To  obtain  values  of 
conversion  efficiency  significantly  greater  than  the  value  0.3  which  Is 
potentially  available  from  Schottky  diode  mixers,  an  Increase  In  I^  with 

a corresponding  sacrifice  In  T Is  necessary.  Note  that  there  Is  no 

in 

fundamental  principle  prohibiting  values  of  r|  greater  than  unity  - the  IF 
power  Is  derived  from  the  bias  supply  rather  than  the  signal. 

From  Fig.  9 It  Is  seen  that  r^^  can  be  substantially  Increased  by 
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detuning  the  circuit  resonance  slightly.  In  this  case  we  should  replace  i 

the  denominator  of  Eq.  (11)  by  (R  + R + X^,  where  X Is  the  external 

o ITx 

reactance  Introduced.  The  Increase  In  r^^  can  offset  the  Increase  In  the 

denominator  due  to  finite  X,  and  yield  Improved  conversion  efficiency. 

However,  the  mixer  noise  temperature  Is  Invariably  degraded  by  such 

detuning,  so  that  an  overall  system  noise  Improvement  will  result  only 

If  the  IF  amplifier  noise  contribution  Is  significant. 

The  values  of  8^  . used  In  preparing  Fig.  10  were  measured  near 
min 

zero  IF  frequency.  Our  nolse-measurlng  system  can  be  used  to  determine 

the  noise  spectrum  for  frequencies  as  large  as  V /lO.  In  contrast  to 

the  current-driven  case,  we  generally  observe  some  frequency  dependence 

of  the  output  noise  for  strongly  coupled  junctions.  For  IF  frequencies 

within  the  coupling  bandwidth,  3^  . rarely  varies  by  more  than  a factor 

min 

of  2 and  most  often  Increases  with  Increasing  frequency.  This  means  that 

the  high  IF  frequencies  which  are  often  used  in  Schottky  diode  heterodyne  i 

systems  should  not  be  necessary  with  Josephson  mixers.  ' 

VIII.  ARRAYS  OF  JUNCTIONS  I 

35 

It  has  been  suggested  that  the  problem  of  Impedance  matching  to 

Q 

the  low  resistance  of  presently  available  thin  film  bridge  junctions 
can  be  alleviated  by  using  series  arrays  of  such  devices.  In  the  current- 
driven  limit.  It  can  be  shown  that  neither  the  conversion  efficiency  nor 
the  noise  performance  of  such  an  array  is  improved  over  what  can  be 
achieved  (In  principle)  with  a single  junction.  The  Impedance  matching 


difficulties,  however,  are  removed  at  both  the  signal  and  the  IF  frequency. 

In  this  respect,  use  of  an  array  is  equivalent  to  impedance  transforming 

a single  junction.  One  exception  to  this  equivalence  is  the  saturation 

power  of  an  array,  which  Increases  linearly  with  the  number  of  junctions. 

It  is  interesting  to  speculate  whether  the  excess  noise  due  to 

hysteresis  effects  in  the  strong  coupling  limit  would  still  exist  in 
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arrays.  It  can  be  argued  that  even  when  the  array  is  collectively 

well  coupled  to  an  rf  source,  the  individual  junctions  are  presented 

with  an  external  Impedance  much  greater  than  their  own  resistance,  and 

so  are  effectively  current  biased. 

To  examine  this  idea,  we  propose  a simple  mechanism  to  explain  the 

instability  of  a bias  point  between  steps  for  a strongly  coupled  single 
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junction:  We  have  found,  in  agreement  with  a previous  study,  that  the 
Impedance  at  the  LO  frequency  of  a current-biased  junction  always  decreases 
with  increasing  bias  voltage.  With  a finite  LX)  source  Impedance,  it  follows 
that  the  LO  current  must  increase  with  increasing  voltage  bias.  However, 
an  increase  of  LO  current  will  further  suppress  the  critical  current. 

The  overall  result  in  the  absence  of  noise  is  a dc  negative  resistance, 
which  prevents  a stable  bias.  It  appears  that  this  mechanism  should 
apply  equally  to  an  array  of  junctions.  Unlike  the  case  of  a single 
junction,  however,  the  thermal  noise  in  the  junctions  is  not  expected 
to  have  a significant  smearing  effect  on  the  overall  I-V  characteristic 
of  an  array.  Hence  the  negative  resistance  will  persitteven  for  fairly 
large  values  of  F.  The  problems  with  strong  coupling  to  an  array  will 
probably  involve  achievement  of  a stable  bias  rather  than  excess  noise. 

The  junctions  in  an  array  must  be  sufficiently  alike  in  their  pro- 
perties that  most  of  them  can  be  biased  near  the  optimum  voltage  = hv^^/Ae. 
If  we  assume  that  the  fabrication  technique  results  in  junctions  with  some 

spread  in  I but  uniform  values  of  RI  , the  variation  in  bias  voltage  <'5V> 
c c 

of  an  individual  junction  from  the  average  <V>  « hV^^/Ae  is  given  by: 
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Here  6l^  is  the  deviation  of  the  critical  current  of  the  junction  from 
the  average,  and  and  are  the  dc  bias  current  and  rf  current, 
common  to  all  the  junctions.  The  term  in  parentheses  is  close  to  unity 
for  optimum  bias  conditions  at  all  values  of  normalized  frequency.  If 
we  compute  from  Eq.  (8)  with  • 100  yA,  T * 4.2  K,  f2  « 0.5,  we  find 
(SV/<V>  w 40  6l/I^.  Hence  only  those  junctions  in  an  array  whose  critical 
currents  are  within  "v*  1%  of  the  average  will  contribute  an  appreciable 
IF  voltage. 


IX.  REAL  JUNCTIONS 


The  theory  developed  in  this  paper  is  relevant  to  real  junctions 
only  to  the  extent  that  they  are  described  by  the  resistively  shunted 
junction  model.  The  earliest  theoretical  justification  for  this  model 
was  supplied  by  Aslamazov  and  Larkin^  for  the  case  of  a three-dimensional 
constriction  joining  two  bulk  superconducting  regions.  If  the  character- 
istic dimension  £ of  the  constriction  satisfies  the  constraint  I « a « ^ 

(Z  is  the  electron  mean  free  path,  C is  the  coherence  length) , the  Ginzberg- 
Landau  (GL)  equations  predict  a sinusoidal  static  current-phase  relation, 
with  the  critical  current  being  related  to  the  normal  resistance  R of  the 
constriction  by  Eq.  (1).  Likharev  and  Yakobson  have  extended  these 
calculations  to  the  case  of  a thin  filament  of  arbitrary  length  connecting 
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two  bulk  pieces  of  superconductor.  They  find  that  as  long  as  the  length 

of  the  bridge  Is  less  than  “V  twice  the  coherence  length,  the  current- 

phase-  relation  does  not  differ  by  more  than  'v<  20%  from  sinusoidal,  and 

the  RI  product  can  acttially  exceed  the  tunnel  Junction  limit  by  'v  30%. 
c 

According  to  the  RSJ  model,  voltage  spikes  with  a width  'v<  Tj  = h/(2eRI^) 
will  occur  under  virtually  any  bias  condition  involving  a dc  voltage 
component.  The  foregoing  analyses  might  not  be  expected  to  be  valid 
unless  the  GL  order  parameter  can  respond  at  least  as  rapidly  as  the 
voltage  across  the  junction  varies.  In  fact,  the  order  parameter  relax- 
ation time  is  predicted  to  be  from  6-12  times  greater  than  Tj  in  those 
instances  where  a time-dependent  GL  equation  can  be  derived  from  the 

38 

microscopic  theory  of  superconductivity.  Recently  Likharev  and  Yakobson 

39 

and  Baratoff  and  Kramer  obtained  solutions  to  the  time-dependent  equations  j 

! 

for  the  case  of  the  thin  filament  joining  two  superconductors.  For  suffi-  ■ 

( 

ciently  short  lengths,  it  is  found  that  the  RSJ  model  is  recovered  with 
corrections  that  amount  to  a small  change  in  the  effective  normal  resis- 
tance and  a cos  ip  modulated  conductance. 

It  is  possible  that  the  dimensional  requirements  for  validity  of  the 
RSJ  model  might  be  met  in  point  contacts,  at  least  if  the  bias  voltage  is 
not  too  large.  However,  with  the  dimensions  that  can  reasonably  be 
expected  to  be  achievable  for  film  type  bridges  (length  comparable  to 
the  coherence  length)  deviations  from  the  RSJ  model  are  apparent  in  the 
calculations.  For  example,  the  voltage  oscillations  are  smaller  in 
amplitude  and  somewhat  skewed.  The  resemblance  is  still  strong  enough 
that  one  might  expect  the  RSJ  model  predictions  of  mixing,  with  an 
"effective",  smaller,  critical  current,  to  be  useful. 


. ■■  r..--  !!:■«  t . . ■_r. 
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The  time-dependent  GL  theory  does  predict  one  very  coimnon  observation  In 
point  contacts  and  bridge  structures:  the  offset  voltage  (or  current) 
that  persists  up  to  very  high  bias  currents.  On  the  other  hand.  It 
cannot  account  for  the  gap  structure  that  is  frequently  seen  in  these 
structures. 

Recently  small  area  tunnel  junctions  with  very  high  current  density 
40 

have  been  reported  to  have  I-V  curves  showing  very  little  hysteresis, 

41 

presumably  because  3^  is  sufficiently  small.  McDonald  £t  al.  have 
reported  the  results  of  calculations  based  on  the  tunneling  formalism  in 

29 

the  limit  of  small  capacitance.  It  has  been  proposed  that  this  model 

should  apply  also  to  high  resistance  point  contacts  whose  contact  area 

is  smaller  than  an  electron  mean  free  path.  The  calculated  I-V  curves 

show  significant  gap  structure,  as  well  as  a much  greater  sensitivity  to 

3^  than  is  characteristic  of  the  RSJ  model.  It  is  plausible  that,  with 

a small  amount  of  noise  rounding,  single-valued  I-V  curves  would  result, 

making  such  junctions  useable  as  mixers.  It  would  be  of  considerable 

interest  to  extend  these  calculations  to  the  case  of  a low  capacitance 

tunnel  junction  biased  with  an  rf  current  source,  allowing  comparison 

with  the  predictions  of  the  RSJ  model. 

Another  feature  of  real  junctions  that  has  not  been  accounted  for 

in  our  analysis  is  the  effect  of  self-heating.  This  has  recently  been 

42 

addressed  by  Tinkham  ^ for  the  case  of  a superconducting  constric- 

tion. He  obtains  a rather  general  expression  for  the  maximum  temperature 
in  the  bridge: 


(29) 
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We  might  guess  that  the  effective  junction  temperature  Is  the  average 

of  T'  and  the  ambient  temperature  T.  Assuming  a bias  point  V = hv/4e 

for  mixer  operation,  it  is  easily  demonstrated  that  quantum  noise  effects 

will  outweigh  the  Increase  in  temperature  implied  by  Eq.  (29). 

A much  more  Important  effect  of  self-heating  is  a suppression  of 

the  critical  current.  Tinkham  et  al.  estimated  this  effect  by  using 

the  time- Independent  GL  equations  with  a spatially  varying  temperature. 

They  found  I '''  exp(-P/P  ),  where  P is  the  power  dissipated  in  the  junc- 
c o 

tlon  and  P^  is  a parameter  that  depends  only  on  materials  parameters  of 
the  superconductor  and  bridge  geometry, 

P C(0)/p  . (30) 

o c 

Here  p is  the  resistivity  of  the  metal.  Nximerically,  it  Is  found  that 

for  most  pure  superconductors  P * 10  yW. 

^ o 

We  can  estimate  the  consequence  of  self-heating  on  I-V  curves  by 
supposing  that  at  finite  voltage  the  RSJ  model  is  obeyed  with  a critical 
current  which  depends  on  the  power  dissipation  in  the  junction.  This 
leads  to  hysteretlc  I-V  curves,  with  the  amount  of  hysteresis  being 
determined  by  6 B RI^/P^.  The  amount  of  hysteresis  Al  between  the  zeroth 
and  first  steps  for  a junction  biased  with  an  rf  current  can  be  calculated 
to  first  order  In  6 using  published  dependences  of  step  heights  for  the 
RSJ  model.^^ If  the  rf  current  suppresses  by  50%  in  the  absence 
of  self-heating,  we  find 


Al  < 0.45  I nfi 
c 


(31) 
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for  0.1  < < 1.0.  For  most  experimental  conditions  a fractional 

hysteresis  AI/I^  * 0.02  would  be  so  thoroughly  masked  by  noise  as  to  be 

unimportant.  At  * 0.3  this  would  require  6 < 0.15.  Using  the  value 

= 10  yW,  and  assuming  an  RI^  product  of  1 mV,  we  find  that  heating 

effects  are  unimportant  if  R > 0.7  ohm.  Empirically  determined  values 

of  for  point  contacts  more  often  lie  in  the  range  1-3  yW,  which  still 

presents  no  serious  limitation  for  practical  mixers.  However,  the  use 

of  alloys  is  certainly  precluded  from  this  point  of  view,  since  Eq.  (30) 

— 3/2 

implies  a dependence  I on  the  electron  mean  free  path. 

CONCLUSION 

We  have  found  that  the  RSJ  model  predictions  of  mixing  performance 

fall  naturally  into  two  regimes:  weak  and  strong  coupling  to  the  signal 

source.  In  the  former  case  the  predictions  take  a particularly  simple 

form.  It  is  found  that  the  mixer  noise  temperature  is  roughly  independent 

of  normalized  frequency  if  0.1<n<1.0,  and  scales  linearly  with  the 

(actual  or  effective)  junction  temperature. 

The  weak  coupling  limit  will  probably  apply  in  most  experiments 

involving  the  high  end  of  the  frequency  range  we  have  considered  if  only 

because  of  practical  limitations.  Our  results  predict  quite  respectable 

performance  nonetheless.  For  instance,  at  300  GHz  7 K)  with 

r = 4,  we  would  have  T 300  K using  a junction  whose  RI  product 
o m c 


exceeds  600  yV.  The  conversion  efficiency  would  be  as  good  as  the  best 
that  could  be  expected  from  a Schottky  mixer. 
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At  lower  frequencies  where  strong  coupling  could  in  principle  be 
achieved,  the  complexity  of  the  observed  behavior  in  our  simulationij 
precludes  definite  predictions.  The  main  overall  conclusion  to  be  drawn 
is  that  for  a fixed  RI^  product  the  normal  resistance  should  be  as  large 
as  possible  - i.e.,  noise  rounding  should  be  pronounced.  This  is  in 
contrast  to  the  weakly  coupl'^d  case,  where  all  that  matters  in  determining 
T^  is  the  ratio  of  junction  and  source  resistances.  The  simulated  results  I 

in  the  strong-coupled  case  are  probably  more  sensitive  to  small  deviations 
from  the  RSJ  model  than  those  for  the  weak-coupled  case.  In  particular, 
the  general  lowering  of  r^^  that  is  implied  in  the  common  observation  of 
"excess  current"  could  result  in  a lower  threshold  in  R for  the  appearance 
of  hysteresis.  This  in  turn  could  reduce  the  excess  noise  in  the  junction, 
resulting  in  a lower  noise  temperature  than  our  simulations  would  predict. 

Finally,  we  emphasize  the  practical  importance  of  saturation  effects 
in  Josephson  mixers.  In  most  applications  the  bandwidth  in  which  the  mixer 
is  coupled  will  have  to  be  restricted.  The  microwave  coupling  circuit 
could,  of  course,  also  function  as  a filter.  Indeed,  an  advantage  of 
using  low  resistance  junctions  is  that  the  impedance  transformation 
required,  if  done  resonantly,  necessarily  implies  a narrow  coupling 
bandwidth.  Only  with  large  arrays  is  the  saturation  level  high  enough 
to  permit  broad  band  coupling. 
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FIGURE  CAPTIONS 


Fig.  1.  The  equivalent  circuit  of  (a)  the  junction,  and  (b)  the  rf  and 
IF  coupling,  that  were  assumed  in  the  simulations.  The  output 
circuit  (c)  of  the  junction  was  found  to  be  valid  for  IF  fre- 
quencies up  to  approximately  1/10  of  the  LO  frequency. 


Fig.  2.  Examples  of  static  I-V  curves  generated  by  a junction  simulator 

for  various  values  of  normalized  frequency,  rf  current,  and  noise 
parameter  F.  The  dependence  of  the  zero-voltage  critical  current 


on  rf  current  shown  at  the  bottom  of  the  figure  is  reproduced  from 
Ref.  15. 


Fig.  3.  The  dependence  of  several  (defined  in  Eq.  (7))  on  normalized 


frequency. 


Fig.  4.  The  dependence  on  normalized  frequency  of  the  parameters  P,  Q, 


and  r^,  used  in  Eq,  (22)  to  obtain  the  differential  resistance. 


Fig.  5.  A comparison  of  the  conversion  efficiencies  that  would  be  expected 
for  various  orders  N of  harmonic  mixing,  for  a range  of  normalized 
signal  frequencies.  It  is  assumed  that  the  LO  frequency  is  near 


V /N.  The  values  chosen  for  r and  F are  plausible  design  goals; 
Slg  b 


the  results  can  be  extended  to  different  values  of  F using  Eq.  (8) 


Fig.  6.  (Dashed  lines)  Simulator  predictions  of  noise  parameter  8*  vs. 


bias  voltage,  for  a junction  with  and  without  LO  current  applied. 


Fig.  6. 
(cont. ) 

Fig.  7. 

Fig.  8. 


Fig.  9. 


Fig.  10. 


Also  shown  are  the  corresponding  current-voltage  curves  (solid 
lines) . A subsequent  recallbratlon  revealed  that  the  noise 
results  shown  are  too  high  by  'Vi  30%. 

The  dependence  of  the  noise  parameters  ^^min^^l 

on  normalized  frequency,  as  deduced  with  the  simulator. 

A typical  sequence  of  I-V  curves  for  a junction  strongly  coupled 
to  the  LO,  using  progressively  larger  values  of  driving  noise  F. 
The  results  are  insensitive  to  the  value  of  the  (unloaded)  Q of 
the  rf  coupling  circuit. 

A demonstration  of  the  effect  on  Junction  I-V  curves  of  detuning 
the  rf  coupling  circuit.  The  LO  level  has  been  adjusted  in  each 
case  to  yield  the  same  suppression  of  the  critical  current. 

A summary  of  simulator  results  for  the  case  of  a strongly  coupled 
junction.  The  mixer  noise  temperature  and  conversion  efficiency 
are  shown  as  functions  of  operating  temperature,  for  the  following 
junction  parameters? 

(A)  1 ■ 110  pA,  R (ohms)  = 0.1  v (GHz) 

c 

(B)  I =•  11  PA,  R (ohms)  = 1.0  V (GHz) 

c 

(C)  I - no  PA,  R (ohms)  - 0.04  V (GHz) 

c 

(D)  I - 11  PA,  R (ohms)  - 0.4  v (GHz) 

c 

Note  that  (A)  and  (B)  correspond  to  ■ 0,19,  while  (C)  and  (D) 
correspond  to  ■ 0,47,  The  unloaded  Q of  the  coupling  circuit 


Fig.  10. 
(cont.) 


was  20  In  each  case.  The  results  are  fairly  sensitive  to 
detuning  of  the  circuit , but  not  to  its  Q when  operated  at  its 


center  frequency.  Thermal  driving  noise  was  assumed  throughout. 


(o)  Junction  Model 


rf  circuit  Junction  I F and  DC  circuit 

(b)  Externol  Coupling  Model 


(c)  Equivalent  Low-Frequency  Circuit 
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